(N 

o 



CO 
CN 



CO 



First principle electronic, structural, elastic, and optical properties of strontium 

titanate 

Chincdu E. Ekuma 1 , Mark Jarrell 1 , Juana Moreno 1 , and Diola Bagayoko 2 

1 Department of Physics and Astronomy & Center for Computation and Technology, 

Louisiana State University, Baton Rouge, Louisiana 70803, USA 

2 Department of Physics, Southern University and A&M College, Baton Rouge, Louisiana 70813, USA 

(Dated: March 1, 2013) 

We report self-consistent ab-initio electronic, structural, elastic, and optical properties of cubic 
SrTiO;? perovskite. Our non-relativistic calculations employed a generalized gradient approximation 
(GGA) potential and the linear combination of atomic orbitals (LCAO) formalism. The distinctive 
feature of our computations stem from solving self-consistently the system of equations describing 
the GGA, using the Bagayoko-Zhao- Williams (BZW) method. Our results are in agreement with 
experimental ones where the later are available. In particular, our theoretical, indirect band gap of 
3.24 eV, at the experimental lattice constant of 3.91 A, is in excellent agreement with experiment. 
Our predicted, equilibrium lattice constant is 3.92 A, with a corresponding indirect band gap of 
3.21 eV and bulk modulus of 183 GPa. 

PACS numbers: 71.15.Mb, 71.20.-b, 71.20.-b, 71.20.Mq, 71.20.Nr 
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I. INTRODUCTION 

Strontium titanate (SrTiOs) is one of the most stud- 
ied oxides of the ABO3 perovskite type structures, due 
to its great technological importance. Many interest- 
ing phenomena such as colossal magnetoresistance, high- 
T c superconductivity, multifcrroicity, and ferroelectric- 
ity are observed in complex oxides. Since most of the 
interesting complex oxides have perovskite structure, 
SrTiOs is an ideal starting point for their study. It 
has been widely used for integration with other oxides 
into heterostructures. Those heterostructres show in- 
teresting properties such as thermoelectricity [1, 2] and 
superconductivity [3, 4]. Many new concepts of mod- 
ern condensed matter and the physics of phase transi- 
tions have been developed while investigating this unique 
material [5-7]. SrTiOs has applications in the fields of fer- 
roelectricity, optoelectronics and macroelectronics. It is 
used as a substrate for the epitaxial growth of high tem- 
perature superconductors. SrTiC>3 exhibits a very large 
dielectric constant. In comparison with SiC>2, SrTiC>3 has 
almost two orders of magnitude higher dielectric constant 
and may as well offer a better replacement for Si02 in 
Si-based nanoelectronic devices (see Wilk et al. [8] [and 
references therein]). SrTiC>3 has found usage in optical 
switches, grain-boundary barrier layer capacitors, cat- 
alytic activators, waveguides, laser frequency doubling, 
high capacity computer memory cells, oxygen gas sen- 
sors, semiconductivity, etc [9-16]. 

During the last few decades, the electronic, structural, 
clastic, and optical properties of SrTiOs (STO), as a 
model of ABO3 perovskite, have been under intensive in- 
vestigation both experimentally[17-23] and theoretically 
[9-13, 24-27]. But, from a theoretical point of view, a 
proper description of its electronic properties is still an 
area of active research. Theoretical computations have 
had difficulty in predicting the correct band gap energy 
and other related electronic properties of SrTiOs from 
first principle. The density functional theory plus addi- 
tional Couloumb interactions (DFT+U) formalism [28- 



31] has had good successes in obtaining correct energy 
bands and gaps of materials, but can only be applied 
to correlated and localized electrons, e.g., 3d or 4/ in 
transition and rare-earth oxides. The hybrid function- 
als (for e.g., Heyd-Scuseria-Ernzcrhof (HSE) hybrid func- 
tional [32-34]) has also been used in attempt to improve 
on the energy bands and band gaps of materials. This 
approach involves a range separation of the exchange 
energy into some fraction of nonlocal Hartree-Fock ex- 
change potential and a fraction of local spin density ap- 
proximation (LSDA) or generalized gradient approxima- 
tion (GGA) exchange potential. Wc should note that 
this range separation is not universal. There is always 
a range separation parameter ui which varies between 
and 00. While it is reasonably clear that there ex- 
ists a value of u that gives the correct gap for a given 
system, this uj is not universal as it is always adjusted 
from one system to another [35, 36]. For example, in 
HSE06 [33, 37], u = O.IIoq 1 (a is the Bohr radius) and 
in Perdew-Burkc-Ernzcrhof (PBEh) global hybrid [38], 
it is 25 % short-range exact exchange and 75 % short- 
range PBE exchange. Even though the HSE functional, 
in most cases, accurately reproduces the optical gap in 
semiconductors, it severely underestimates the gap in in- 
sulators [36, 39] and its band width in metallic systems 
is generally too large [36, 39-41]. The Engel and Vosko 
[42] (EV) GGA and the Tran and Blaha [43] modified 
Bcckc- Johnson (TB-mBJ) have also provided some im- 
provements to the band gap of materials. For TB-mBJ, 
while the band gaps arc considerably improved, the effec- 
tive masses are severely underestimated [41]. In the case 
of the EV potential, the equilibrium lattice constants are 
far too large as compared to experiment and, as such, 
leads to an unsatisfactory total energy [44, 45]. 

The theoretical underestimations of band gaps and 
other energy eigenvalues have been ascribed to the inad- 
equacies of density functional potentials for the descrip- 
tion of ground state electronic properties of semiconduc- 
tors [18, 19, 24]. Also, other methods [46-48] that en- 
tirely go beyond the density functional theory (DFT) do 
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FIG. 1: (Color online). The iso-surface cubic unit cell of 
SrTi03 with the iso-lines, at lattice parameter of 3.91 A. 
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FIG. 2: Calculated, band structure of c-SrTi03 at the ex- 
perimental lattice constant of 3.91 A as obtained with the 
optimal basis set using PW - GGA. The horizontal, dashed 
line indicates the position of the Fermi energy (Ejr) (-0.10872 
eV) which has been set equal to zero. 



not obtain the correct band gap values of most semicon- 
ductors without adjustment or fitting parameters [49, 50]. 
This unsatisfactory situation is a key motivation for our 
work. 

Stoichiometric STO has an experimental, indirect band 
gap of 3.20 - 3.25 eV at room temperature [17, 19, 20]. 
Theoretical calculations using several techniques have led 
to band gaps of SrTiOs in the ranges 1.71 to 2.2 eV for 
LDA and GGA [6, 10, 12], 1.87 to 3.63 eV for Hybrid 
DFT [5, 9, 10], and value as high as 11.97 eV for the 
Hatree-Fock (HF) method [10]. 

In this paper, we present a simple, yet robust, and ab- 
initio method, based on self consistent solutions of the 
pertinent system of equations [51-54], that correctly pre- 
dicts band gap values and related electronic properties of 
SrTiOs rigorously, from first principle, within the LCAO- 
GGA formalism. We also compute the structural, elastic 
and optical properties of SrTiOs. 

The rest of this article is organized as follows. After 
this introduction in section I, the computational meth- 
ods and details are given in section II. The results of our 
self-consistent calculations are presented and discussed in 
section III. We then summarize and conclude in section 
IV. 



II. METHODS AND COMPUTATIONAL 
DETAILS 

In the ground state, STO has the simple cubic (0\ - 
Pm3m) perovskitc structure [55], with Sr atom sitting at 
the origin, Ti at the body center and three oxygen atoms 
at the three face centers [56] (see Fig. 1). We used the 
room temperature experimental lattice constant of 3.91 



A [55, 56]. 

Our ab initio, self consistent, nonrclativistic calcula- 
tions employed a linear combination of atomic orbitals 
(LCAO) formalism and a generalized gradient approxi- 
mation (GGA) potential. One may argue that relativistic 
effects are important for the description of SrTiOs. As 
was noted by Marques et al. [57], relativistic effects are 
only important for the description of the high-K dielectric 
band structure. Their calculated, relativistic and non- 
rclativistic band structure for SrTiOs, up to an energy of 
5 eV for the valence and conduction bands, respectively, 
are almost identical. Consequently, we expect a negligi- 
ble relativistic correction for the band gap of SrTiOs. 

The distinctive feature of our calculations, the use of 
the Bagayoko, Zhao, and Williams (BZW) method, has 
been extensively described in the literature [51-54, 58- 
60]. This method has been shown to lead to accurate 
ground state properties of many semiconductors: c-InN 
[60], w-InN [54], w-CdS [53], c-CdS [61], rutile-Ti0 2 [59], 
AlAs [62], GaN, Si, C, Ru02 [51], and carbon nanotubes 
[63]. 

Instead of assuming that a single trial basis set will 
yield the correct ground state charge density of the 
solid, the BZW method entails a minimum of three self- 
consistent calculations with basis sets of different sizes 
and generally with different polarization functions, i.e., 
p, d, and f functions. The correct ground state is the 
one where all the occupied energies are at their minima. 
In practice, up to seven self-consistent calculations have 
been performed for some materials e.g, wurtzite ZnO [64]. 
The computations begin with a relatively small basis set 
that should not be smaller than the minimal basis set. 
The latter is defined as one that is just large enough to 
account for all the electrons in the atomic or ionic species 
present in the solid. The preliminary, self-consistent cal- 
culations of the electronic properties of the species pro- 



vide the wave functions that serve as input in the solid 
state calculations. 

The first, self consistent calculation for the solid is per- 
formed with this small basis set (Calculation I) that is 
subsequently augmented with one orbital for the next 
self-consistent calculation (Calculation II). Depending on 
the angular symmetry of the added orbital, the size of 
the new basis set is larger than that of the previous one 
by 2, 6, 10, or 14 for s, p, d, and f functions, respec- 
tively. The occupied energies from calculations I and II 
are compared graphically and numerically. For the first 
two calculations, we found these occupied energies to be 
different for all the solids we have studied to date [65, 66], 
including SrTiC>3. The basis set for calculation II is then 
augmented in order to carry out self-consistent calcula- 
tion III. Again, the occupied energies from calculations 
II and III are compared. This process of augmenting the 
basis set and of performing self-consistent calculations is 
continued until the occupied energies from a calculation, 
say TV, are found to be identical to their corresponding 
ones from calculation (N+l), within computational un- 
certainties that are less than 50 meV. This perfect super- 
position of occupied energies from two consecutive calcu- 
lations identifies the basis set for Calculation N as the 
optimal one, i.e, the smallest basis set that yields the 
lowest, occupied energies of the system. The attaine- 
ment of this minima signifies that this basis set is verifi- 
ably complete for the description of the ground state of 
the system. Larger basis sets that include the optimal 
one do not lower any occupied energies from their values 
obtained with the optimal basis set. 

As explained elsewhere [51, 53], these larger basis sets 
do not lead to any changes in the ground state charge 
density or the Hamiltonian. However, larger basis sets 
that include the optimal one often lead to a lowering of 
some unoccupied energies. This lowering of some un- 
occupied energies cannot be ascribed to a physical in- 
teraction included in the Hamiltonian. Up to the opti- 
mal basis set, changes in the basis sets lead to changes 
in the charge density, the potential, and the Hamilto- 
nian. Hence, changes in occupied and unoccupied en- 
ergies, for self-consistent calculations with basis sets up 
to the optimal one, can be ascribed to a physical inter- 
action embedded in the Kohn-Sham Hamiltonian. The 
system of equations in DFT totally determines changes 
in the occupied states. It also determines, at least in 
part, low-laying unoccupied states that are interacting 
with the occupied ones, up to the optimal basis set. For 
example, in wurtzite InN, the calculated dielectric func- 
tions agree with their corresponding experimental ones 
up to 5.5-6.0 eV [67]. Given that only direct transitions 
were taken into account in our dielectric functions cal- 
culations, this agreement indicates that the low-laying 
unoccupied bands were correctly determined. For larger 
basis sets that include the optimal one, the extra lower- 
ing of some unoccupied energies is a direct consequence 
of the Raylcigh theorem [53, 58, 59, 65]. This theorem 
states that when an eigenvalue equation is solved with 
basis sets I and II, with set II larger than I and including 
I entirely, then the eigenvalues obtained with set II arc 



lower than or equal to their corresponding ones obtained 
with basis set I [53, 59]. 

The above process entails iterations for the equation 
giving the ground state charge density, with the itera- 
tions for the Kohn-Sham equation carried out for each 
choice of the basis set. Given that iterations for the 
Kohn-Sham (KS) equation involves the charge density 
(CD) equation, one could conclude that a single trial ba- 
sis set calculation solves both equations self consistently. 
A problem with this view stems from the fact that any 
two such calculations, with different trial basis sets, lead 
to different, converged (i.e., self consistent) eigenvalues of 
the KS equation. The fundamental theorem of algebra 
suffices to guarantee that the two sets will be different if 
the basis sets have different numbers of basis functions 
utilized in the expansions. The question then arises as 
to which of the two sets of eigenvalues of the KS equa- 
tion provides the physical description of the system un- 
der study. To answer such a question definitively and 
from first principle, the BZW method follows the pro- 
cess described above. Upon reaching the optimal basis 
set, not only the charge density, the potential, and the 
Hamiltonian no longer change (i.e., they have converged), 
but also the resulting, self-consistent eigenvalues of the 
KS equation have reached their respective minima, for 
the occupied states. In our understanding, to solve the 
system of equation sclf-consistcntly means obtaining con- 
verged eigenvalues (attainable with most arbitrary basis 
set) but also occupied eigenvalues that have reached their 
respective minima (attainable with BZW method). 

In the above sense, the BZW method solves self- 
consistently not just the Kohn-Sham equation, but also 
the equation giving the ground state charge density in 
terms of the wave functions of the occupied states. Fur- 
ther, in his Nobel lecture [68], Kohn noted the "den- 
sity optimal" feature of the wave functions from correctly 
performed DFT calculations while those for the Hartree 
Fock approach are "total-energy optimal." Without a 
constrained search for the converged ground state, it is 
quite difficult to infer the basis set that yields the correct 
ground state charge density [69, 70]. This point becomes 
clearer by noting that the reorganization of the cloud 
of valence electrons is drastically different for atomic or 
ionic species as compared to molecules or solids. For in- 
stance, single trial basis set calculations cannot make up 
for any deficiency in the angular symmetry of the func- 
tions, irrespective of the degree of convergence of the it- 
erations of the Kohn-Sham equation. By correct ground 
state charge density, we mean the charge density that 
leads to the minima of all the occupied energies. 

In this work, we utilize the electronic structure pack- 
age from the Ames Laboratory of the US Department of 
Energy (DOE), Ames, Iowa [71]. We employ the gener- 
alized gradient approximation (GGA) potential given by 
Perdew and Wang [72] . We utilize sets of even tempered 
Gaussian functions with exponents from 0.12 to 10 5 to 
form the atomic wavefunctions. There are 15, 15, and 13 
s, p, and d orbitals, respectively, for Sr, while 17, 17, and 
15 s, p, and d orbitals, respectively, are used for both Ti 
and O. The charge fitting error using the Gaussian func- 



tions in the atomic calculation is about 10 . Since the 
deep core states are fully occupied and are inactive chem- 
ically, the charge densities of the deep core states were 
kept the same as in the free atom. However, the core 
states of low binding energy were still allowed to fully 
relax, along with the valence states, in the self consistent 
calculations. The orbitals employed in the self-consistent 
calculations are between paranthesis for Sr (3d 4s Ap Ad 
5 s), Ti (3p 3d 4s) and O (2s 2p 3s), including some that 
are unoccupied in the free atoms or ions. These unoc- 
cupied orbitals arc included in the self-consistent LCAO 
calculation to allow the restructuring of the electronic 
cloud, including possible polarization, in the crystal en- 
vironment. 

The Brillouin zone (BZ) integration for the charge den- 
sity in the self consistent procedure is based on 56 special 
k points in the irreducible Brillouin zone (IBZ). The com- 
putational error for the valence charge is 5.3 x 10 -5 eV 
per valence electron. The self consistent potential con- 
verged to a difference of 10 -5 after several tenths of it- 
erations. The energy eigenvalues and eigenfunctions are 
then solved at 161 special k points in the IBZ for the 



band structure. A total of 152 weighted k points, cho- 
sen along the high symmetry lines in the IBZ of SrTiC>3, 
are used to solve for the energy eigenvalues from which 
the electron density of states (DOS) are calculated using 
the linear analytical tetrahedron method [73]. The par- 
tial density of states (pDOS) and the effective charge at 
each atomic sites are evaluated using the Mulliken charge 
analysis procedure [74]. We also calculated, the equilib- 
rium lattice constant a , the bulk modulus (B ), the as- 
sociated total energy and the electron and hole effective 
masses in different directions. 

In calculating the lattice constant, we utilize a least 
square fit of our data to the Murnaghan's equation of 
state [75, 76]. The lattice constant for the minimum total 
energy is the equilibrium one. The bulk modulus {B ) is 
calculated at the equilibrium lattice constant. 

The dielectric function e(cj) = e\{oS) + i£2(w) can be 
calculated once the electronic wave functions and ener- 
gies are known. The imaginary part of the dielectric 
function £2(w), from the direct interband transitions, is 
calculated using the Kubo-Greenwood formula [77]: 
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where tvjj is the photon energy, P = —ihV is the momen- 
tum operator, Q is the volume of the unit cell, ipknif) and 
ipkiif) are the initial and final states, respectively, ff-i is 
the Fermi distribution function for the i t h states, and e^i 
is the energy of the electron in the i t h state. The real 
part of the dielectric function, £i(w), is obtained from 
the well-known Kramers-Kronig (KK) relation, 



ei(w) = l + -Af 
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du', 
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where M indicates the principal value of the integral [78] . 
The real part of the optical conductivity, Re[a(u>)], fol- 
lows from above as 

Re[a(w)]) = ^£ 2 M (3) 



III. RESULTS AND DISCUSSION 

The results of the electronic structure computations 
are given in Figs. 2 to 4. Figure 5 depicts the calcu- 
lated total energy of STO, while Fig. 6 shows the optical 
spectra obtained using the optimal basis set from the 
electronic structure computations. Figures 1 and 4 have 
been drawn using the xcrysden [79] . We discuss the elec- 
tronic structure, low laying conduction bands, and the 
effective mass in III A. The DOS is discussed elaborately 
in IIIB. The structural properties are presented in IIIC, 
while III D deals with the calculated, optical properties. 



A. The Electronic Structure, Band Gap, 
Low-energy Conduction Band and Effective Mass 

The electronic structure of the valence and the low en- 
ergy conduction states determine the band gap and other 
important properties of materials. Table I shows that our 
ab initio, first principle method yielded an indirect band 
gap of 3.24 eV at the experimental lattice constant of 
3.91 A (see Fig. 2) and an indirect band gap of 3.21 eV 
at the calculated equilibrium lattice constant of 3.92 A. 
This table contains several, previous results from calcu- 
lations using LDA, GGA or hybrid potentials. Table II 
contains the calculated energies at some high symmetry 
points in the BZ. These energies are provided for future 
comparison with experimental and theoretical findings. 
Our calculated band structure (see Fig. 2) resembles that 
of the parent Ti02 system. Fig. 2 also shows that the top 
of the valence band is at the L point. 

The effective mass is one of the main factors deter- 
mining the transport properties, the Seebeck coefficient, 
and electrical conductivity of materials. The calculated 
electron effective masses at the bottom of the conduction 
band along the V - L, V - X, and V - K directions, respec- 
tively, are 0.68 - 0.81, 0.44 - 0.59, and 0.51 - 0.66 while the 
calculated hole effective masses at the top of the valence 
band, along the Y - L, T - X, and T - K directions, respec- 
tively, are 0.64 - 0.83, 1.22 - 1.27, and 0.96 - 1.02 (all in 
units of the electron mass). The observed anisotropy and 
the ranges of effective masses confirmed the earlier obser- 
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FIG. 3: (Color online), (a) Calculated, density of states (DOS) of c-SrTiOs at the experimental lattice constant of 3.91 A 
as obtained with the optimal basis set using PW - GGA. (b) Calculated, partial density of states (pDOS) of c-SrTi03 at the 
exnerimenta.l lattice constant of 3.91 A as nhta.ined with the optimal basis set using PW - GGA. The vertical, dashed line 

been set equal to zero. 










FIG. 4: (Color online). The contour plot of the electron 
charge density of c-SrTi03 as obtained with the optimal basis 
set using PW - GGA. The center atom is Sr with Ti atoms at 
the corners and O atoms in the middle of the size of the par- 
allelogram around Sr. A n(r) is the variation of the electron 
charge density as a function of distance away from an atomic 
site. A logarithmic scale has been used. 



vations of Mattheiss and co-workers [80, 81]. Our calcu- 
lated effective masses are in excellent agreement with the- 
detailed effective mass values as reported by Mattheiss 
[80] [and references therein] and the relativistic computa- 
tional results of Marques et al. [57]. 

There is a significant 2p - T^hd hybridization in the 
valence bands. As shown in Fig. 2, the valence bands 
of SrTiC>3 can be divided into three distinct groups: the 
upper, intermediate, and lower groups of valence bands 
occupy the energy ranges of to -5.80 cV, -14.2 to - 
14.78, and -16.62 to -17.80 cV, respectively. The upper 
VB bands are made up of nine bands with a bandwidth 
of 5.80 cV, in agreement with X-ray photocmission spec- 
troscopy (XPS) values of 5 - 6 cV [871 and the GGA 




FIG. 5: (Color online) Calculated, total energy per unit cell 
as a function of the lattice constant of c-SrTi03, as obtained 
with the optimal basis set using PW - GGA. The calculated 
equilibrium lattice constant is 3.92 A. 



results of Jiangni ct al. [12]. They are formed by the hy- 
bridization between O 2p and Ti 3d, with very little con- 
tribution from the Ti 3p and Sr 4p (see Fig. 3(b)). Two 
of the bands at the T point are triply degenerate: rig 
(-2.75 eV) and 2 ri 5 (-0.25 eV), while the third band T 25 
(-1.23 eV) is non-degenerate (see Table II). Immediately 
below the upper VB, a group of triply degenerate bands 
emanating from the hybridization between Sr 4p, O 2s 
and O 2p, with little dispersion at 3 Fi 5 . This group is 
located at -14.78 eV. The lowest lying VB bands are semi 
core like bands formed mainly due to the hybridization 
between O 2s, and Sr 4p with very little bonding coming 
from Ti 4s and Ti 3p. They are located at ^ (-17.80 
eV) and 2 Ti (-16.62 eV). Our calculated position of Sr 
Ap between - 14.78 to - 17.80 cV is in agreement with the 
XPS measurement of Battye ct al. [88] which placed it 
at 16.50 eV below the Fermi surface (Ep). Also, Board 
ct al. [89] and Battye et al. [88] measured the average 



TABLE I: Comparison of our calculated band gap values with 
other theoretical and experimental ones, for c-SrTi03. Our 
calculations show that c-SrTi03 has an indirect band gap 
from L to V points. Exactly the same band gap value is 
found from L to X points. All the band gaps are indirect 
unless otherwise indicated with (D). 





Computational Method 


E s (eV) 


GGA 


GGA 


- BZW (Present work) 


3.24 




(with equilibrium lattice constant) 


3.21 






PP - PWGA 


1.97 a 






PP - PBE 


1.99 a 






FP - LAPW 


1.80 6 






PP 


1.60 c 






FLAPW 


1.78 d 






TB - LMTO 


1.40 (D) e 






FB - LMTO 


2.20 (D) 7 


LDA 




PP 


2.04 a 






PP- PW 


1.79 m 






LMTO - ASA 


1.80 9 






PP 


1.71 ft 






PP- PW 


1.79 m 






OLCAO 


1.45" 


HF 


PP 


11.97" 


Hybrid DFT 




B3PW - LCAO 


3.63* 






PP - BLYP 


1.94 a 






PP - B3LYP 


3.57 a 






PP - P3PW 


3.63 a 






LCGO - B3PW 


3.70 J 


SA 


FP - LAPW 


1.87 - 3.25 fc 


Experiment 


NA 


3.10 - 3.25' 



a Ref. 


[10 


6 Ref. 


[61. 


c Ref. 


12 


d Ref. 


[57 


e Ref. 


[82 


•fRef. 


[S3 


»Ref. 


[84 


'•Ref. 


[49 


"Ref. 


[91. 


•*Ref. 


[51. 


fc Ref. 


[13 


'Ref. 


[56 


'"Ref. 


[85 


"Ref. 


[86 



position of the O 2s to be about 17 eV below Ep. Our 
calculated value of -17.85 eV in Fig. 3(a) is close to the 
experimental one. In particular, our result does not un- 
derestimate this core state position as is often the case 
in GGA calculations. 

The conduction bands (CB), immediately above the 
Fermi level (low energy conduction band) , are dominated 
by threefold degenerate Ti 3d £ 2 g orbitals which hybridize 
with O 2p and O 2s. The two- fold degenerate Ti 3d e g 
states have some hybridization from all other orbitals ex- 
cept Ti (4s and 3p), Sr (4p and 3d). The energy eigen- 
value in the lowest conduction bands, at the X point, are 



practically the same as that at the T point, resulting in 
the observed, minimal dispersion in the conduction band 
between F and X points. This feature is apparent from 
our energies at the high symmetry points in Table II. At 
the T point, energies associated with the lowest-laying 
conduction bands are: ^25' (3.24 eV), 1 r 12 (4.84 eV), 
2 r 12 (6.10 eV), 3 T 1 (6.75 eV) and 2 T 2b , (8.92 eV). The 
calculated, low energy conduction bands in Fig. 2 are 
quite different from those of previous studies. Figure 2 
shows that the lowest conduction bands are degenerate 
at the r and X points along the [100] and equivalent 
directions. The electronic structure in Fig. 2 was calcu- 
lated using the experimental lattice constant. Wc further 
examined whether or not the position of the shallow mini- 
mum in the lowest conduction band depends on the value 
of the lattice constant by using several values of the lat- 
tice constant around the experimental one. Even though, 
the band gap value changed from 3.26 to 3.17 eV, there 
was no appreciable change in the depth of the shallow 
minimum in the lowest conduction band. We recall that 
the gap is 3.21 eV for our calculated, equilibrium lattice 
parameter. 



TABLE II: Eigenvalues (eV), along high symmetry points, for 
c-SrTi03, as obtained with the experimental lattice constant 
of 3.91 A. The Fermi energy of - 0.12188 eV is set to zero 
in the table. The energy eigenvalues at F and X points are 
found to be almost identical. 



L 


r 


X 


K 


-32.06 


-32.03 


-32.03 


-32.07 


-32.06 


-32.03 


-32.03 


-32.04 


-32.06 


-32.03 


-32.03 


-32.03 


-16.19 


-16.62 


-16.62 


-16.65 


-16.19 


-16.62 


-16.62 


-16.49 


-14.96 


-14.67 


-14.67 


-14.61 


-14.96 


-14.67 


-14.67 


-14.47 


-14.96 


-14.67 


-14.67 


-14.33 


-4.94 


-2.75 


-2.75 


-3.92 


-4.32 


-2.75 


-2.75 


-3.43 


-4.32 


-2.75 


-2.75 


-3.09 


-3.58 


-1.23 


-1.23 


-1.68 





-0.25 


-0.25 


-1.33 





-0.25 


-0.25 


-1.04 





-0.25 


-0.25 


-0.73 


5.22 


3.24 


3.24 


4.06 


5.22 


3.24 


3.24 


4.12 


5.22 


3.24 


3.24 


4.33 


8.65 


4.82 


4.82 


5.38 


10.92 


10.28 


10.28 


13.42 


10.92 


10.28 


10.28 


14.47 



TABLE III: Experimental and theoretical lattice constants a (in 
modulus (in GPa). 



for c-SrTiC>3 along with the calculated values of the bulk 





Computational Method 


a (A) 


B (GPa) 


GGA 


BZW - LCAO (Present work) 


3.92 


183.45 




PP - PWGA 


3.95 (3.93) 


167 (195)° 




PP - PBE 


3.94 (3.93) 


169 (195)" 




PW 


3.95 (3.88) 


167 (194) 6 




PBE 


3.95 


167 (194) c 




11 


3.91 (3.82) 


210.21 (252.92) d 


LDA 


PP 


3.86 


214 (215)° 




LAPW 


3.86 


204 (176) c 




FLAPW 


3.95 


167 c 




•>•> 


3.93 (3.87) 


207.28 (227.63) d 




PP- PW 


3.87 


194 e 




OLCAO 


3.93 


163 / 


HF 


PP 


3.92 (3.93) 


219 (211)° 




PP 


3.98 (3.92) 


208.85 (206.68) d 


Hybrid DFT 


PP - BLYP 


3.98 


164° 




PP - B3LYP 


3.94 


177 (187)" 




PP - P3PW 


3.90 (3.91) 


177 (186)° 


Experiment 


NA 


3.89 - 3.92 9 


174 - 183 ft 



°Ref. 


10]. 


6 Ref. 


90]. 


c Ref. 


11]. 


d Ref. 


91]. 


e Ref. 


85]. 


•*Ref. 


86]. 


9Ref. 


55, 56, 92, 93] 


h Rcf. 


56, 92, 93]. 



B. Densities of States, Electron Distribution and 
Chemical Bonding 



Figures 3(a) and 3(b) exhibit the total (DOS) and 
related partial (pDOS) densities of states, respectively. 
Figure 4 shows the contour plot of the distribution of 
the electron charge density of SrTi03. As can be seen 
from Fig. 4, the electron density of SrTi03, away from 
the atomic sites, does not have a spherically symmetric 
distribution. Further, the bonding between Ti and O is 
covalent, due to Ti-3d and 0-2p hybridization, unlike in 
the case of Sr and O that is ionic. The bond length of 
Sr - O is 2.76 A, with a minimum charge density of ~ 
0.19 e/A 3 , while Ti - O bond length is 1.95 A, with a 
charge density of ~ 0.63 e/A 3 . The experimental bond 
lengths for Sr - O and Ti - O are 2.76 and 1.96 A, re- 
spectively, [94, 95] with corresponding charge densities 
of 0.2 and 0.67 - 0.90 e/A 3 , in that order [94-97]. The 
charge density distribution around the O atom with re- 
spect to the horizontal Ti - O - Ti line is elongated in the 
direction along the Ti - O covalent bond in agreement 
with room temperature experimental charge density dis- 
tribution of SrTiC>3 reported by Ikeda et al. [96]. This 
anisotropic charge density distribution is ascribed to the 
rotational mode of the Ti - 06 octahedron by experiment 
[94, 95, 981. 



From our calculated DOS (see Fig. 3(a)), it can be in- 
ferred that the onset of absorption is quite sharp and it 
starts at about 3.24 eV. It exhibits a fine structure at 3.6 
eV, with a shoulder at 4.50 eV. This picture is consistent 
with the experimental results of Cohen and Blunt [99] 
and the theoretical findings of Perkins and Winter [100] 
of a relatively sharp absorption edge in the optical mea- 
surement of SrTiOa. In the calculated DOS of the low 
laying conduction bands, sharp peaks appear at 4.70 
eV, 5.72 eV, and 7.05 eV. Relatively broad peaks are 
found at 8.33 eV, 10.97 eV and 12.75 eV. Our com- 
puted peaks are in basic agreement with experimental 
findings of Cardona [101] and Braun et al. [102]. For 
the valence bands DOS, we calculated peaks at -0.20 
eV, -0.72 eV, -1.14 eV, -1.83 eV, -2.86 eV, -4.06 eV, - 
4.50 eV, -14.30 eV, -4.61 eV, and -17.53 eV. The peaks 
in the valence band DOS are all sharp. Our calculated 
electronic structure is in agreement with scanning trans- 
mission electron microscopy, vacuum ultraviolet spec- 
troscopy, and spectroscopic ellipsometry measurement of 
Van Bcnthcm et al. [103]. 

Our calculated band gap value of 3.24 eV, from L to 
r, is practical the same as the experimental one. In gen- 
eral, other theoretical calculations obtained values that 
are up to 1.1 eV smaller. The source of the small band 
gap values was believed to be due to the pushing up of 
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FIG. 6: (Color online), (a) Calculated, dispersive part, £i(a;), of the dielectric function of c-SrTi03 as obtained with the optimal 
basis set using PW - GGA. (b) Calculated, absorptive part, £2(1^), of the dielectric function of c-SrTi03 as obtained with the 
optimal basis set using PW - GGA. (c) Calculated, optical conductivity, ct(oj), of c-SrTi03 as obtained with the optimal basis 
set using PW - GGA. As per the inserts (see Figs. 6(b) and 6(c)), the band absorption edge is at 3.24 eV. In all cases, spectra 
have been calculated without any broadening. 



the top of the valence band dominated by Ti 3d and O 
2p states to higher energies [22]. According to our find- 
ings, it rather appears to be the extra lowering of the 
conduction bands that produces GGA (or LDA) band 
gaps that are more than 1.1 cV smaller than the experi- 
mental values, if LCAO type computations do not search 
for and utilize an optimal basis set. Such a basis set is 
verifiably converged for the description of occupied states 
[53, 60, 61]. We recall that Kohn and Sham [104], in their 
original paper, explicitly stated the need to solve self con- 
sistently the pertinent system of equations defining LDA. 
The BZW method, as explained above, rigorously solves 
the system of equations in the sense explained above. 
Single trial basis set calculations also involve both the 
KS and charge density equations. The major difference 
resides in the fact that these calculations do not gener- 
ally entail changes in the basis functions beyond those of 
the expansion coefficients. 



C. Structural Properties 



The total energy versus the lattice constant data are 
shown in Fig. 5. The data fit well to the Murnaghan 
equation of state (EOS). The calculated equilibrium lat- 
tice constant is 3.92 A and the bulk modulus, B , is 
183.45 GPa. 

The experimently reported lattice constants are in the 
range 3.89 to 3.92 A [55, 56, 86, 92] and the bulk modulus 
lays in the range 174 to 183 GPa [56, 86, 92]. In Table III, 
we show our calculated equilibrium lattice constant and 
bulk modulus in comparison with experimental and other 
theoretical results. Both our calculated lattice constant 
and bulk modulus agree well with corresponding, exper- 
imental ones, respectively. 



D. Optical Properties 

The plot of the dispersive (ei(w)) and absorptive 
(e 2 (w)) dielectric functions are shown in Figs. 6(a) 
and 6(b), respectively, while the optical conductivity 
(tr(o;)) profile is in Fig. 6(c). All reported spectra have 
been calculated without any broadening and, may have 
more features than experimental ones. Our calculated, 
dielectric spectra are in good agreement with the exper- 
imental measurements [93, 101]. The calculated optical 
spectra only included the direct inter band transitions. 
The fundamental absorption edge E , which is also a 
measure of the optical gap, was found to be 3.24 eV from 
the calculation, as per the insert of Fig. 6(b). Our com- 
puted direct gap of 3.43 eV is in agreement with the 
experimental value of 3.40 cV[93, 101]. Our calculated 
£i(w = 0) at zero frequency equals 4.75 (cf. Fig. 6(a)). 
It compares well with the experimental value of 4.92 mea- 
sured by Braun et al. [102]. Our dielectric spectra re- 
semble that of BaTi03 of Bagayoko et al. [105]. This 
observation also holds for the data of Cardona [101] and 
Baurerlc et al. [106]. Both the experimental and our cal- 
culated results show that the direct optical gap is larger 
than the smallest indirect band gap. 

Figure 6(c) shows the optical conductivity o-(uj) of 
SrTiOa. As per the insert, it also shows that the funda- 
mental absorption edge starts at 3.24 eV. The positions 
of the peaks (without any rigid shift) are in agreement 
with experimental findings. 



IV. SUMMARY AND CONCLUSION 

We have performed first principle, ab initio calcula- 
tions of the electronic, structural, elastic, and optical 
properties of bulk SrTiOa in the cubic phase using GGA 
potential and the BZW method. Our calculated results, 
without any adjustment or corrections, show good agree- 
ment with experimental data. 



The agreement of our calculated band gaps (3.21 and 
3.24 eV) and electron effective masses with correspond- 
ing, experimental values is significant. Some calculations 
with adjustable parameters can lead to the correct band 
gap; but they generally do not yield the correct curva- 
ture of the conduction band around its minimum-as given 
by the electron effective masses. Similarly, the agreement 
between the peaks in the calculated density of states with 
corresponding, experimental ones denote the correct de- 
scription of the relative location of the bands. This result 
is confirmed by our reproduction of the measured features 
of the dielectric functions, the imaginary part of which 
was obtained using only direct transitions between occu- 
pied and unoccupied bands. Our calculated equilibrium 



lattice constant (3.92 A) and bulk modulus (183.45 GPa) 
also agree with corresponding, experimental findings. 
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